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Abstract 

\ir  Force  requirements,  such  as  the  knowledge  of  the  upper  atmosphere  en¬ 
vironment  of  vehicles  and  the  knowledge  of  the  propagation  characteristics  of  radio 
and  radar  signals,  require  the  solutions  of  the  equations  of  motion  of  fluid  dynam¬ 
ics  and  of  electromagnetic  theory  which  are  often  very  complicated.  This  report 
presents  a  new  mathematical  approach  to  the  obtaining  of  such  solutio  is.  The 
vector  field  is  represented  in  such  a  form  that  new  techniques  may  be  used  to  find 
the  appropriate  solutions.  Come  problems  of  fluid  dynamics  and  electromagnetic 
theory  are  solved  as  an  illustration  of  the  new  approach.  In  later  reports,  new 
techniques  will  be  used  in  other  problems  related  to  Air  Force  needs. 

In  this  report,  eigenfunctions  of  the  curl  operator  are  introduced.  The  ex¬ 
pansion  of  vector  fields  in  terms  of  these  eigenfunctions  leads  to  a  decomposition 
of  such  fields  into  three  modes,  one  of  which  corresponds  to  an  irrotational  vector 
field,  and  two  of  which  correspond  to  rotational  circularly  vector  fields  of  opposite 
signs  of  polarization.  Under  a  rotation  of  coordinates,  the  three  modes  which  art 
introduced  in  this  fashion  remain  invariant.  Hence  the  Helmholtz  decomposition 
of  vector  fields  has  been  introduced  in  an  irreducible,  rotationally  invariant  form. 

These  expansions  enable  one  to  handle  the  curl  and  divergence  operators 
simply.  As  illustrations  of  the  use  of  the  curl  eigenfunctions,  four  problems  are 
solved. 
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Eigenfunctions  of  the 
Curl  Operator,  Rotation  ally  Invariant 
Helmholtz  Theorem,  and  Applications  to 
Electromagnetic  Theory  and  Fluid  Dynamics 


1.  INTRODUCTION 

In  problems  which  involve  vector  tields,  the  curl  and  divergence  operators 
play  a  central  role.  For  example,  in  electromagnetic  theory.  Maxwell's  equation 
are  expressed  in  terms  of  these  operators,  while  in  fluid  dynamics  it  is  convenient 
to  decompose  the  velocity  vector  field  into  two  parts:  one  for  which  the  curl  Is 
zero,  and  the  other  for  which  the  divergence  is  zero.  In  this  report,  eigenfunctions 
of  the  curl  operator  are  introduced  which  considerably  simplify  problems  involving 
vector  fields,  particularly  when  these  fields  are  defined  over  the  entire  coordinate 
space.  Vector  fields  in  terms  of  these  eigenfunctions  will  be  expanded,  and  will 
show  that  each  vector  of  a  vector  field  can  be  expressed  as  a  sum  of  three  vectors: 
one  is  an  irrotational  vector,  and  the  other  two  are  rotational  vectors  which  are 
circularly  polarized  and  have  opposite  signs  of  polarization.  Furthermore,  this 
decomposition  is  rotationally  invariant  because  in  a  frame  of  reference  obtained 
from  the  original  one  by  a  rotation,  the  three  vectors,  into  which  the  rotated  vector 
is  split,  each  maps  separately  into  the  corresponding  vector  into  which  the  original 
vector  was  decomposed.  That  is,  the  irrotational  vector  goes  Into  the  irrotational 
vector,  while  each  of  the  two  rotational  circularly-polarized  vectors  go  into  ‘he 
corresponding  circularly-polarized  vector.  This  allows  improvement  upon  the 
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original  Helmholtz  decomposition  of  a  vector  into  a  rotational  and  irrotational 
part.  Vector  and  scalar  potentials  will  also  be  introduced  in  a  rotationally  in¬ 
variant  manner  euch  that  the  irrotational  mode  is  obtained  from  a  scalar  potential, 
and  each  of  the  circularly-polarized  rotational  modes  is  obtained  as  the  cur)  of  a 
corresponding  vector  potential.  It  will  also  be  shown  how  to  construct  these 
vector  and  scalar  potentials  from  the  vector. 

In  addition  to  giving  the  expansion  of  the  vectors  in  terms  of  cartesian 
variables  for  the  arguments  of  the  vector  functions,  expansions  will  be  given  in 
terms  of  spherical  coordinates.  In  this  case,  the  modes  are  conveniently  ex¬ 
pressed  in  terms  oi  the  vector  spherical  harmonics.  This  expansion  will  obviously 
be  useful  for  cases  of  spherical  symmetry. 

rhe  initial  value  problem  for  Maxwell's  equations  will  be  solved  when  the 
sources  and  currents  are  given  functions  of  time  and  apace,  and  will  show  how  the 
use  of  the  curl  eigenfunctions  yields  ill  a  very  simple  manner  the  uncoupling  of  the 
longitudinal  and  transverse  modes.  The  velocity  vector  fields  of  fluid  mechanics 
can  also  be  interpreted  in  terms  of  these  modes  and  show  that  the  irrotational 
inodes  of  the  vector  field  correspond  to  irrotational  motion  of  the  fluid,  while  the 
circularly-polarized  components  correspond  to  the  rotational  motion  of  the  fluid. 
The  vorticity  can  be  related  nimplv  to  these  circularly-polarized  components. 

Such  concepts  will  be  used  to  obtain  families  of  new  exact  solutions  of  the  in¬ 
compressible  Navier-Stokes  equations  and  solutions  in  the  theory  of  linearized 
compressible  viscous  flow.  Exact  solutions  will  also  be  obtained  for  vertical 
Shear  motion  in  an  ocean  on  the  rotating  earth. 

The  eigenfunctions  of  the  curl  and  their  properties  were  obtained  by  consider¬ 
ing  helicity  representations  of  the  rotation  group.  The  general  properties  of  such 
representations  have  been  discussed  in  great  detail  by  Moses  (1970).  The  method 
by  which  the  eigenfunctions  were  obtained  is  given  in  Appendix  A.  Having  obtained 
the  eigenfunctions,  however,  it  is  a  straightforward  matter  to  verify  their  proper¬ 
ties  by  direct  computation,  and  this  verification  will  be  indicated  in  the  body  of 
this  repor*. 


a.  THE  qa  vectors 

Let  the  variable  X  take  on  the  three  values  X  «  1,  0,  -1,  and  Q,be  any  unit 
vector.  Then  for  every  fa  r02‘  hj)  the  three  vectors  ”  l  <  H.). 

QjxOl  >.  Q^fg.)]  we  h. 


20(a>  *  -a.  ■ 


Qx(n)  *  -x( 2 


i"  1/2 


rJ1(i?1+Un2) 

r+n. 


rJ2(r).+  lAr)  ) 

-*■  — i^r'~-iX-ni+iA,|2 


a) 


for  A.  *  ±  1  . 

The  following  properties  of  are  easily  verified: 


gx^a*-  2„<a>  a*  ■ 


£Qix(2.)cVa>  *  4ij  - 

(2) 

a.  ■  g  *  o  for  x  ■  ±  i , 

a  x  (a*  *  *^Sx  <a*  ■ 

rhe  above  equations  show  mat  the  unit  vectors  span  the  three -dimensional 
vector  space  for  each  vector  jj_  • 

The  vectors  ^  were  introduced  previously  (.n  a  somewhat  different  notation) 
as  a  means  of  reducing  the  electromagnetic  fields  and  the  electromagnetic  vector 
potential  to  the  Irreducible  representations  of  the  Poincare*  group  (Moses.  19fifia). 
This  report  extends  the  applicability  of  these  vectors. 

The  following  properties  of  Qx  are  useful  when  it  is  desired  to  impose  reality 
conditions  upon  the  vector  fields: 

Sx  '  <a>  ’  -  for  X  «  i  1, 

So ;  (a.>  “  9o(a}  ■  <3) 


Also 


Let  be  any  vector.  Then  g^(£>  is  defined  by 

Sx<a>52x<a>- 


(4) 


(5) 
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* 


whore 


n  * 


(5a) 


Then 


■  -f  ■ 


yA(p)  »  ->.(2  ) 


•1/2 


P1(P1'fUP'2) 

p(pvp3) 


-  1, 


p.^p^Up^) 

p(p+p3) 


(0) 


for  X  *  ±  1  . 


It  will  now  be  convenient  to  militate  how  the  vectors  Q  A  behave  under  rotations 
of  the  frame  of  reference.  As  is  well  known,  every  rotation  of  a  frame  of  refer¬ 
ence  can  be  described  by  means  of  a  vector  d  1  (d.,  ll.t,  d.()  where  the  direction  of 
0  gives  the  axis  of  rotation  and  d  «  |  <)  |  is  the  angle  of  rotation  (0  1  rt  t).  If  a 
vector  x^  lias  the  three  components  x;  in  one  frame  of  reference,  in  a  rotated  frame 
it  will  have  the  components  x.1  where 


V  *  £  • 


(7) 


and  K^(£)  are  the  matrix  elements  of  the  rotation  mati  ix  H  (il)  . 


Hij(d)  .  6ijCos  d  ♦  0td 


+  ^  ”"ijk,,k  d 


sin  0 


(8) 


where  e^  ts  the  usual  anti-symmetric  icnsur: 


*jik  "'ikj  ’  -eijk'  *12  3  *  1  • 


(8a) 


It  is  convenient  to  define  the  vector  x  '  as  being  the  vector  whose  components 
In  the  unrotated  frame  are  given  by  x.'  which  are  the  components  of  in  the 
rotaw.  frame.  Then  Eq.  (7)  may  be  written  as 


x  '  ■  H(0)x  . 


(7a  > 


5 


Let  be  a  vector  function  of  jm  Then  under  a  rotation  parametrized  by 

0  .  the  components  of  this  vector  in  ‘be  rotated  frame  will  be  components  of  the 
vector  Vi, 1  ( Jt, ’ )  given  by 

Jt'(x)  ■  R(£)ii[H(-£)jL]  •  (7b) 

In  a  similar  way,  one  can  define  vector  functions  of  the  vector  p  and  give 
their  transformation  properties  in  tfms  of  the  vector  which  describes  a  rota¬ 
tion  of  coordinates.  The  vectors  Qx  (^)  may  be  regarded  as  being  vector  functions 
of  .  In  a  certain  sense  they  are  "eigenfunctions"  of  rotations-  They  satisfy  the 
relation 


R(0)gx  (R(-2.)^J  =  exp{2iX4>(?,n)}Qx  (£>  ■ 


(9) 


where 

P 


(9a  > 


and  ®  (ft ,  n  )  is  the  principal  branch  of 
(&'  0g)tan(O/2) 

Un  ®  *  fJf'l+rjg)  +  <£  -  a)73  UnTo'/'a)  ■  1 

The  proof  of  Eqs.  (9)  and  (Ob)  will  be  given  in  Section  5,  where  additional 
properties  of  the  vectors  will  be  given. 


3.  EIGKNK1  NOTIONS  01  THU  U  hl.  OPKIUTOH.  KAPWSiON 
THEOREMS,  \NU  INN  \RI\NT  IIKI.MIIOI.)  Z  THEOREM 


Let  the  vectors  _x^(x_|  be  defined  by 

*A  <*  Ig.)  *  <2*)‘3/VC-'  •  00) 


Then  from  the  properties  of  the  Fourier  transformation  and  the  completeness 
and  orthogonality  properties  of  the  vectors  §  x  in  Kq.  (2),  the  vectors  X  X(x  |j^) 
satisfy  the  following  orthogonality  and  completeness  relations: 


* 


/dJL  X\ ■<£  lc_)  >  =  6>m6<£  LL’^  • 

£fdPxi\'(llZ)x  =  V—’*' 


(11) 


where  x  denotes  the  i'th  component  of  the  vector  X  A 
From  Eq.  (2)  it  is  seen  that 


CvX^(x|p)  *  P*  Os.||u). 


V  -  x  x  (x|)^)  *  0  for  *  *  i  1  .  ( 12 ) 

V  •*  j  (5.1k.)  ~  -  ip  (2  3,,2r’£.^- 


fhe  first  part  of  Eq.  (12)  shows  that  the  vectors  X  A  are  eigenfunctions  of 
the  curl  operator  with  eigenvalue  Ap.  The  variable  X  itself  may  be  considered  to 
be  the  eigenvalue  of  the  operator  (-£2)  v  when  this  operator  is  properly 

interpreted. 

The  completeness  and  orthogonality  relations  of  Eq.  (11)  enable  us  to  expand 
any  vector  function  u  (x^),  subject  to  very  general  restrictions,  as  follows: 


u(x)  «  £  J  dp  xA(x|p)gA(p) 


*  lu^lx) 


(13) 


where 


Hi*** 


*x(^|p)8x<£)  ■ 


(13a) 


Eqs.  (13)  and  (13a)  constitute  our  sharpening  of  the  Helmholtz  decomposition 
theorem,  for  from  Eq.  (2) 


V-uA(x^)“0.  forA»±l, 
V  X  u  -  (x  >  «  0  . 

-X.  0 


(14) 


One  can  also  introduce  vector  and  scalar  potentials.  For  X  «  ±  I,  the  vector 
potentials  A^(x)  and  the  scalar  potential  V(x_)  are  introduced.  First  note  that  the 
coefficients  g  ^ (jo )  of  the  expansion  Eq.  (13)  can  be  btalned  fro~«  u(x)  as  follows: 

6x<E.>  •/  x  ^  *  (  x|  j>)  *  i^(  £)  dx_ .  (13b) 

Then 

£x<*>*  x/{t)  *>(~|P)gX(P)  • 

V(x  )  -  *  gQ  (p)  .  (15) 

From  Eq.  (12)  it  follows  that 
3£_)  *  V  x  ^^(x  ),  for  X  ■  ±  1 

u  0(x)  *  VV(x)  ■  (16) 

Thus  the  original  Helmholtz  theorem  has  been  sharpened  in  two  ways.  First, 
two  irrotational  vec  jrs  have  been  introduced  in  the  decomposition  of  a  general 
vector,  each  of  which  is  the  curl  of  its  own  vector  potential.  Second,  a  procedure 
lias  been  given  for  obtaining  the  vector  and  scalar  potentials  when  the  vector  u  (x ) 
is  given.  Now  it  will  be  shown  how  our  formulation  is  rotationally  invariant. 

Under  the  rotation  of  axes  described  by  the  vector  6,  the  components  of  the 
vectors  (x)  are  given  in  the  new  frame  by  the  components  of  u^’(x  ),  where 

u  X’(x)  -  H(S)u  x[R(-0)xl  ,  (17) 

See  Eqs.  (7a)  and  (7b). 

On  using  Eq.  (9),  note  that  g^’  (p),  which  replaces  gj^(g_)  in  the  expansion 
Eqs.  (13)  and  (13a),  is  related  to  g^(g_)  by 

gx’(|3>  *  exp{2iX$(£,  rj)  }gxf  R(-  0>£j  .  (18) 


where  the  rnction  <t(0_,  q)  is  given  by  Eq.  (9b). 
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Thus,  exct  pt  for  a  phase,  the  decomposition  given  by  Eqs.  (13),  (13a),  and 
(13b)  Is  invariant.  Vectors  belonging  to  values  of  X  do  not  mix,  unlike  the  situa¬ 
tion  for  the  cartesian  components. 

The  scalar  potential  and  the  vector  potentials,  which  in  the  rotated  frame 
are  denoted  by  V!(x_)  and  A'(x_)  are  also  uncoupled: 

V'(x)  *  V(  R(-  £)x] 

Ax'(x)  =  R{0)  Ax{R(-6)x]  .  (19] 

Vector  and  scalar  potentials  do  not  mix.  Furthermore,  there  are  no  am¬ 
biguities  in  gauge  in  the  present  invariant  set-up. 

Now  it  will  be  shown  how  the  requirement  that  u(x)  be  real  affects  the  ex¬ 
pansion  Eq.  (13)  and  puts  restrictions  on  the  functions  gx(^).  From  Eq.  (4)  it  is 
seen  that  a  necessary  and  sufficient  condition  that  the  vector  u  (x^>  be  real  is 


«A<P>  *  * 


/prixP2\ 

V  Pl+iXp2  / 


gx*<-p) 


When  the  vector  u  (x )  is  real,  it  will  be  convenient  to  replace  the  expansions 
Eqs,  (13),  (13a)  and  (15)  by 

W  -t|e / dE.ix<*  k)«x<p>  +£ /dp  *x*<*  k>e**(p) 

u  X<x)  >  (13c 


where 


«-x(x)  "T  j/dE.  ix<*lE.)8x<P> +/dp  *X*<xIR.>«X  *<P>  • 

Also 

*>(r'  ■  i>'il  t'sx'e1  */(?)  1  x * ‘5. 1 11  *x*<p >  • 

vn.1  (p,  -/(f),-*': 5o.,p,  . 
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».  application," 

This  section  gives  some  applications  of  the  use  of  the  curl  eigenfunctions, 
and  shows  how  the  modes  that  have  been  introduced  uncouple  in  many  cases. 

4.1  Application  to  Maxwell's  Equations 

4. 1. 1  THE  INITIAL  VALUE  PROBLEM  IN  THE  USUAL  FORM. 

In  terms  of  Gaussian  units,  Maxwell's  Equations  are: 


4*  1  (x.t) 

VxH(x,t)x  4c^(x,t)+i-  ~at 

(21) 

,  dH(x,t) 

V  XE(x,t)  =  -  jr  -~gr~ —  - 

(22) 

V  •  E  (x  ,  t)  «  4lrP(x,*  t)  , 

(23) 

Z  ’  H  <x  ,  t)  -  0  . 

(24) 

The  current  density  j^(x_,  t)  and  the  charge  density  p(x  ,  t)  are  taken  as  given 
functions  of  x  and  t,  and  it  is  assumed  that  the  dielectric  constant  e  *  1  and  the 
magnetic  permeability  n  =  1. 

Then  the  general  initial  >  °ilue  problem  posed  by  these  equations  will  be  solved. 
Actually,  in  terms  of  the  present  notation  it  will  be  shown  how  the  modes  uncouple. 
In  paragraph  4.  1.  2,  a  more  convenient,  though  less  conventional,  notation  will  be 
introduced,  and  then  the  solution  of  the  initial  value  problem  will  be  carried  out  to 
the  end. 

Expansion  Eq.  (13)  for  the  vectors  E  and  H  may  be  written  as 
H(x, t)  *  P  x  | 


E  (x ,  t) 


*  k>fx<e.-t)  * 


(25) 


where  the  functions  h^tp,  t)  and  f^(p,t)  are  regarded  as  unknown  functions  of 
their  arguments. 

Also,  the  known  current  density]  is  expanded  as 

jjx.t)  ■  £ /dP  i.x(*  IS.)kA(P’° ’ 


(26) 
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* 


where  the  functions  kx(j^,  t)  art  known  from  Eq.  (13b). 

The  scalar  p  is  expanded  in  terms  of  an  ordinary  Fourier  transformation  as 
follows: 


P<x,t) 


e*S. 


i.r(p.  t)  . 


(27) 


where  the  function  r(p,  t)  is  also  known. 

•N* 

Substitute  Eqs.  (25)  and  (27)  into  Eqs.  (23)  and  (24).  On  using  Eq.  (12),  one 
obtains 


-i(2»)  P  e *£.*£•  hQ(p,  t)  «  0  , 

-i(2r)  3/2/djipe^if0(p,t)»(2^'‘i/dp  eHL’ *  r(p  ,  t)  .  (28) 


Thus  the  longitudinal  components  of  the  electromagnetic  field  can  be  solved 
immediately: 


hQ  (p  .  t)  *  0  . 

V5.- l)  “irr'e.' 0 


(29) 


Similarly,  on  using  Eq.  (25)  and  Eq.  (26)  in  Eqs.  (21)  and  (22).  and  using  the 
fact  that  the  vectors  X^(x_  |  j i)  are  linearly  independent. 


-Sr - 4’ko 


(30) 


and 


1  3hx(p.t) 

Xpfx(p  .  t)  ■  -  -  —gf1 —  .  *  -  ±  1  » 


(31) 


4,  ,  3f>(p.t) 

Xphx(t.  t)  -  —  kx<p  .  t)  +  -  — af .  *-±l 


(32) 


Note  that  the  modes  are  completely  uncoupled.  This  result  is  in  sharp  contrast 
to  the  situation  in  which  one  uses  simply  the  Fourier  transform  of  the  cartesian 
components  of  the  vectors,  where  considerable  coupling  occurs. 


n 


From  the  second  part  of  Eqs.  (29)  and  from  Eq.  (30),  one  obtains  a  necessary 
condition  for  the  solubility  of  Maxwell's  Equations: 


3r(p,t) 

~jt—  ’  *Pk0%.‘ t}  • 


(33) 


This  equation  is  simply  the  equation  of  continuity 
3 p(x ,  t) 

—fif — +  *  ^(x.tJ-O  (33a) 


when  it  is  written  in  terms  of  modes  through  the  substitution  of  p  and  j  by 
Eqs.  (26)  and  (27).  It  is  interesting  to  note  that  only  the  lamellar  part  of  the 
current  density  ^is  determined  by  the  c  large  density  p. 

Since  Eq.  (29)  gives  hQ  (p ,  t)  and  fQ(p  ,  t),  one  need  only  find  h^(p ,  t)  and 
f^(p  ,  t)  for  X  ■  ±  l.  Equations  (31)  and  (32)  are  a  pair  of  ordinary  differential 
equations  which  are  of  first  order  and  which  are  coupled.  They  are  easily  solved 
using  standard,  well-known  techniques.  However,  instead  of  carrying  out  the 
solution  in  the  present  notation.  Maxwell's  Equations  will  be  written  in  terms  of 
another  notation  which  will  give  the  results  more  directly.  The  objective  in  this 
section  was  to  show  the  effectiveness  of  the  use  of  the  curl  eigenfunctions  without 
the  use  of  notational  simplifications.  Furthermore,  it  is  believed  that  the  tech¬ 
niques  used  in  this  section  are  capable  of  being  extended  in  various  directions. 

We  may  go  further  into  the  matter  in  later  reports. 

4.  1.  2  THE  SOLUTION  OF  THE  INITIAL  VALUE  PROBLEM  IN  TERMS  OF 
AN  ALTERNATIVE  NOTATION 


Define  the  complex  vector  ♦  (x  ,  t)  by 
+  (x,  t)  =  E(x,  t)  -  iH(x,t). 


(34) 


In  terms  of'  ♦,  Maxwell's  Equations,  Eqs.  (21)  through  (24)  become 
„  i  3£(x,  t)  4)r. 

!*♦&«- -5- -5t r-  -  irV*-"' 


(35) 


^  ‘  ♦  (x,  t)  =  4jt  p(x,  t)  . 


(36) 


The  expression  ♦  (x,  t)  is  expanded  as  follows: 


(37) 
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Upon  substitution  into  Eqs.  (35)  and  (36)  along  with  Eqs.  (26)  and  (27),  one  obtains 


go(p' 0  *  71 


i  3g0(p.t)  .  , 


(38) 


3gA(p  ,  t) 

- g== - \cp\gx(p  ,  t)  =  -4)rkx(p  ,  t)  ,  for  k  =  ±  1  . 

fhe  first  part  of  Eqs.  (38)  gives  the  irrotationai  or  longitudinal  part  of  the 
electromagnetic  field  immediately.  For  the  second  part  of  Eq.  (38)  to  be  valid, 
the  consistency  condition  Eq.  (30)  is  required  as  in  the  earlier  treatment.  The 
last  part  of  Eq.  (38)  is  an  ordinary  differential  equation  which  gives  the  transverse 
field.  The  general  solution  is 

gA(p,t)  -  gA(p)elcXpt-  4*elcXpt  f  e'lcXpt'kA<^,t,)dt’  ,  (39) 

•'0 

where  gA(p )  is  a  constant  of  integration  and  is  given  by 

gx(p  )  *  •  (40> 

The  first  term  on  the  right  of  Eq.  (39)  represents  the  radiation  field,  and  the 
second  term  represents  the  field  due  to  the  currents.  To  find  the  significance  of 
the  variable  X  in  the  electromagnetic  case,  let  us  consider  the  pure  radiation 
field.  Accordingly,  set  k^(p  ,  t)  =  0  in  Eq.  (39),  and  take  gA(p)  to  be  given  by 

gA(^)  *  6X  16(p1>«(p2)  6(p3-k)  ,  (41) 

where  k  >  0  . 

The  electric  field  is  given  by 

E  j(x,  t)  »  (  ,»)'2/2(2)”  1^2cos(kct+kz), 

E2(x,t)  ■  -{2t)~2/2<2)~ */2sin(kct+kz), 

E3<x.t)  -  0  . 


(42) 
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The  electromagnetic  wave  is  clearly  a  wave  moving  in  the  negative  z-direction 
and  is  circularly  polarized  in  the  direction  of  propagation.  Let  us  now  take 

gx(p>  *  _16(p1)6(p2)6(p3-k)  .  (43) 

The  electric  field  is  given  by 

E.(x,t)  -  -(2^)*3/2(2>'1^2cos(kct-kz). 

E2(x(t)  "  (2>r)‘2^(2)  ^2sin(kct-kz)  (44) 

E3(x  ,  t)  ■  0  . 

This  electromagnetic  wave  propagates  along  the  positive  z-axis  with  a  circular 
polarization  opposite  to  the  direction  of  propagation.  Hence,  in  our  representation, 
the  radiation  field  consists  of  a  superposition  of  circularly-polarized  waves 
travelling  in  various  directions.  Of  course,  by  picking  the  function  gx(p  ) 
appropriately,  plane  polarized  waves  can  also  be  represented. 

4.  1.  3  GREEN'S  FUNCTIONS,  AND  AN  INVERSE  PROBLEM 

By  using  Eq.  (13b),  one  can  express  g-.(p  )  in  terms  of  <Hx,0).  Similarly, 

'*  -V  "V 

kx(p,  t)  can  be  expressed  in  terms  of  the  current  density  j  (x,  t>.  Thi-n,  from  the 
substitution  of  Eq.  (39)  into  Eq.  (37),  it  is  seen  that  ♦  (x,t)  is  an  expression  in¬ 
volving  Green's  functions  in  x,t  space.  For  the  sake  of  brevity,  this  expression 
will  not  be  given  in  this  report. 

It  is  also  possible  to  solve  an  inverse  problem.  Assume  that  the  current 
density  jjx  ,  t)  is  identically  zero  for  t  <  0,  is  "switched  on"  during  the  interval 
0  <  t  <  T,  and  is  identically  zero  again  for  t  >  T.  During  the  times  that  the  cur¬ 
rent  density  if  "off",  the  transverse  part  of  the  electromagnetic  field  will  be  a 
radiation  field.  One  can  use  Eq.  (39)  to  solve  the  direct  problem,  in  which  one 
prescribes  the  radiation  field  before  the  current  is  switched  on  and  calculates  the 
radiation  field  after  the  current  is  switched  off.  The  current  density  is  assumed 
to  be  known.  In  the  inverse  problem,  one  prescribes  the  radiation  field  before 
the  current  density  is  switched  on  and  the  radiation  field  after  the  current  is 
switched  off.  One  asks  for  the  current  density,  during  the  time  that  it  is  "on", 
which  will  take  the  initial  radiation  field  into  the  final  radiation  field.  This  prob¬ 
lem  has  been  shown  to  have  non-unique  solutions  (Moses,  1958)  using  a  different 
form  for  the  electromagnetic  field.  The  present  notation  seems  more  useful  for 
solving  the  inverse  problem.  In  particular,  it  appears  easier  to  impose  "physical" 
conditions  -  such  as  the  requirement  that  the  current  density  be  expressible  in 
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terms  of  multipoles  contained  within  a  given  volume  -  upon  the  inverse  problem 
to  yield  a  meaningful  unique  solution.  Perhaps  we  shall  discuss  this  problem 
in  later  reports. 


1.2  \ppli  ration  lo  fluid  Dynamic. 

4.2.1  IRROTATIONAL  AND  ROTATIONAL  VELOCITY  MODES, 

VORT1CITY  MODES,  INTERPRE  TATION  OF  VOR  TICITY  MODES 
AS  SHEAR  VELOCI  TY  MODES,  AND  STRATIFIED  VORTICITY 
MODES 


I.et  us  now  regard  the  vector  u(x)  which  is  decomposed  into  the  curl  modes 
of  Eqs.  (13)  or  (13c)  as  being  the  velocity  vector  of  a  fluid,  using  the  Eulerian 
description.  The  mode  u^  <x  )  is  clearly  the  irrotationai  mode  and  can  be  ob¬ 
tained  from  a  scalar  potential  V(x)  as  in  Eq.  (lfi).  For  the  present,  only  the 
rotational  modes  u^(x  )  for  which  ^  =  ±  1  will  be  discussed.  From  Eq.  (13c)  it  is 
seen  that  the  most  general  real  rotational  velocity  field  u(x)  can  be  written,  on 
absorbing  some  factors  in  gx(p>  . 


u(x)  «  5- 


P  yde,Sxl2)t,l~~Sx<p)+  P  fd£%\  <ll> 


(45) 


where 


Thus  the  general  rotational  velocity  vector  u(x)  is  a  superposition  of  the 
velocities  uA(x|d.  p."  )  where 


[3,  p,a  )  »  gx(a)eip<H‘  ~a>  +  QA*(rl)e'lp(2.'  i+a). 


(46) 


In  Eq.  (46),  1]  is  any  unit  vector,  p  any  leal  number,  and  a  any  other  real  number. 

The  velocity  uA  (x  |  n,  p,  0)  will  be  called  a  "vorticity  mode"  of  helicity  X,  direc¬ 
tion  n.  wave  length  (2x/p),  and  phase  a.  Note  that  in  a  single  vorticity  mode,  the 
phase  can  be  made  zero  by  changing  the  origin  of  coordinates  appropriately.  To 
find  the  velocity  field  determined  by  a  single  vorticity  mode,  let  us  take  the  phase 
a  ■  0,  and  assume  that  the  unit  vector  r(  points  in  the  positive  z  direction;  that  is, 
ij  »  (0.  0,  1).  Then 

uXl(x!n,  p,  0)  »  <2)1^2X  cos  pz  , 

ux2 (2.1a- p-o)  *  #ln  p*- 


“*3 


|n  ~  A\  -  A 

F'  'Jf  “ 


(47) 


Thus  the  velocity  vector  u^(x  |rj.pa*)  is  transverse  to  the  direction  n.  is 
constant  in  magnitude,  and  rotates  clockwise  as  one  moves  in  the  positive  r) 
direction  of  A  »  l  and  rotates  counter  clockwise  if  A»-l.  Thus,  if  an  axis  is 
assigned  to  the  direction  n  and  if  in  each  plane  orthogonal  to  this  axis  the  velocity, 
vector  is  drawn  from  the  axis,  the  locus  of  the  tip  of  the  vector  will  be  a  negative 
helix  if  A  ■  l  as  one  moves  in  the  positive  direction  on  the  axis,  and  a  positive 
hi  lix  if  A  *  -l.  The  pitch  of  the  helix  equals  the  wave  length  of  the  vorticity  mode, 
namely  'Jn/ p. 

It  is  clear  that  the  vorticity  mode  co.  esponds  to  a  kind  of  shear  velocity 
field.  The  use  of  the  curl  eigenfunctions  leads  to  the  decomposition  of  rotational 
flows  into  shear  flows.  The  notion  of  "little  vortices",  which  one  might  expect  to 
be  useful  in  describing  fluid  motion  with  vorticity,  plays  no  role. 

It  will  be  useful  to  introduce  the  superposition  of  vorticity  modes  which  have 
the  same  direction  o  .  Such  a  superposition  will  be  called  "a  stratified  vorticity 
mode".  Designating  such  a  mode  by  u  (x  |  q), 

£  ^A(H.>  J  dp  eip<IL  V  gx(p>  ,  (48) 

"  x«±l  ~  ■*’ 

where  p  is  a  single  real  variable  and  g^(p)  is  a  complex  function  of  p.  The  vector 
q  is  the  direction  of  the  stratified  vorticity  mode.  One  sees  that  u^(x_|  n  )  as  a 
function  of  x  depends  only  on  x  •  >)  .  It  is  for  this  reason  that  the  .emi  "stratified" 

--V. 

is  used  to  describe  the  superposition.  These  modes  satisfy  the  condition 

1  u  J  u  (x  j  q)  =  0  .  (49) 

As  will  be  seen,  a  consequence  of  liq.  (49)  is  that  one  of  the  non-linear  terms 
in  the  equations  of  fluid  motion  vanishes  and  that  one  is  able  lo  get  some  exact 
solutions  of  fluid  flow. 

From  the  first  part  of  l.q.  (3),  a  necessary  an  i  sufficient  condition  for 
u  (x  |  q)  to  be  real  is 


* 


if. 

On  using  Eq.  (50)  in  Kq.  (48) 

u(x|q)  “  f  rip  olp^L' 

+  Q.W  8p  cip*2.'i*g  :-(p)  ■  (51) 

"**■’  ^  ^  -00 

Let  us  now  define  fx(p)  for  p  >0  by 

>  X<P>  =>  g^P).  f.j  (p)*  -g^'P)  ■  (52> 


Then  on  using  the  first  part  of  Eq.  (3),  Eq.  (51)  becomes 


u(x|n)  -  £  i1  | £>  * 

X*tl 


(53) 


where 


*  2x<l>  Jp  dP  elp(l-i.)fx(p) 

+  Wx*<n>  /"  dP  c-Wl'VfiHp). 


(53a) 


Clearly,  the  vector  u^OcJq)  is  a  superposition  of  vorticity  modes  u  X(x|r^,  p,a) 
belonging  to  the  same  value  of  X  and  n  . 

4.2.2  EXACT  SOLUTIONS  OF  THE  INCOMPRESSIBLE  NAVIER -STOKES 
EQUATION 

4.  2.  2.  1  General  Stratified  Vorticity  Motion 

The  Navier-Stokes  equations  for  an  incompressible  fluid  are 


dv(x  ,  t) 


...  ■  (  v  (x  ,  th  V  I  v  (x  .  t)  +  —  V  P(x  ,  t)  -  i/  v  2v  (x  ,  t)  *  0 

V  ■  v  (x  ,  t)  »  0  . 


(54) 


In  Eq.  (54),  v  <^_,  t)  is  the  velocity  vector,  P(x  ,  t)  is  the  pressure,  v  is  the 
kinematic  viscosity  p/p,  and  p  is  the  density  taken  to  be  constant,  of  course.  The 
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objective  is  to  find  solutions  of  this  equation  which  correspond  to  a  stratified 
vorticity  mode.  That  this  is  possible  is  due  to  the  fact  that  the  non-linear  term 
of  the  first  part  of  Eq.  (54)  vanishes;  see  Eq.  (49),  Furthermore,  it  will  be  seen 
that  the  pressure  P  uncouples  from  the  stratified  vorticity  mode  which  simplifies 
the  problem  even  more.  One  hopes  that  one  can  get  some  insight  into  the  nature 
of  more  general  viscous  flows  and  even,  perhaps,  a  deeper  understanding  of 
turbulence.  In  a  later  report,  wc  shall  give  other  exact  solutions  of  rotational, 
viscous  flow  which,  in  a  sense,  is  complementary  to  the  flow  discussed  in  the 
present  section. 

Actually,  a  somewhat  more  general  problem  than  the  Uine-development  of  a 
stratified  vorticity  mode  will  be  solved.  It  will  be  assumed  that  the  velocity 
v  (x  ,  t)  is  the  sum  of  a  velocity  V  (t),  which  depends  on  time  but  not  on  position, 
and  a  stratified  vorticity  mode.  That  is,  the  stratified  mode  "rides"  a  velocity 
which  varies  with  time  but  whose  space  variation  is  slow  compared  with  the  space 
variation  (that  is.  pitch)  of  the  helices  in  the  stratified  mode. 

Thus 

v  (x  ,  t)  =  V(t)  +  u  (X  ,  0  ,  (55) 


where 


u  (x  ,  t) 


£  2 


dp  e^  ^g.fp,  t>. 


(55a) 


Clearly,  v  (x_,  t)  satisfies  the  second  part  of  Eq.  (54).  Hence  one  need  only  be 
concerned  with  the  first  of  these  equations. 

Let  us  write 

r  :  5V  (t) 

P(x,  t)  «  /  dp  e  Ci.r(p.  t)  -p  - - x  .  (56) 

It  should  be  noted  that,  while  in  Eq.  (55a)  the  intc-gration  is  over  a  smgl<  variable 
denoted  by  p,  the  integration  in  Eq.  (56*  is  over  the  three-dimensional  space  whose 
variables  are  collectively  denoted  by  . 

Using  Eq.  (49)  with  it  (x  ,  t)  for  u  (x  Jr^),  the  nonlinear  term  in  tile  first  part  of 
Eq.  (54)  disappears.  On  substituting  Eqs.  (55),  (55a)  and  (56)  into  the  first  part 
of  Eq.  (54) 


* 


IB 


L  Qx<n)  f^dpc^V 
x«±i 


ag>(p. ') 

-—t - +  ‘-P  h-  v  (t)Jgx(p,  t) 


+  rP  gA(P'  l) 


'  J rip  el£.  i.  pr(p,  t)  . 


(57) 


-ik  ■  x 


Multiply  through  i  y  t  ~  --  and  integrate  with  respect  to  x  : 


A=±l' 


PP  g\<P*t) 


agx(p,t) 

- ^ - +  ipl  >r  V(t)]gx(p.  t) 


=  -  ikr(k.t)  . 


(53) 


t.et  us  dot  through  by  Js  and  use  ky  Q  x^2.^0S."P2.)  =  P  H.'  Q  x^  ^  "PI )  ~  0  • 
Then 


k2r(k,  t)  -  0  (59) 

for  which  the  general  solution  is 

r (k  )  =  A(t)fi(k  )  ,  (BO) 

-w 

where  A (t)  is  on  arbitrary  function  of  t,  as  can  be  shown  from  the  theory  of  dis¬ 
tributions.  Thus  from  Eq.  (56) 

3V(t) 

P(x,  t)  =  A (t >  -  - - x  .  (61) 


Except  for  the  derivative  of  the  space-independent  wind,  the  pressure  un¬ 
couples  from  the  velocity  field. 

Now  substitute  Eq.  (60)  in  Eq.  (57).  The  right  hand  side  of  Eq.  (57)  vanishes. 
On  letting  ?ad_  x.  taking  the  Fourier  transform  with  respect  to  f  ,  and  using 
the  linear  independence  of  the  vectors  gx(r;),  one  obtains  the  following  uncoupled 
ordinary  differential  equations  for  gx(p,  t): 

9gx(P-t)  (  2  I 

+  j  ipl  17  .  V  (t)  J  +  p  r  |  gA  (p,  t)  =  0  , 


3t 


(62) 
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for  which  the  general  solution  ia 


g\<P< t)  “  g^(p)exp 


-ip  f 

•'n 


df  |n  V(t’)j 


exp  [ -vp  tj 


(63) 


where  g^(p)  is  the  constant  of  integration:  g^(p)  =  g>(p,  0)  . 

For  u^(x  ,  t)  to  be  real,  it  ia  required  that 

g.j(p)  *  -gj'M-p)  •  (64 > 

Now  define  f^(p)  for  p  >  0  by 

f X(P>  1  gj(p>.  f.^p)  =  -gj*(-p)  ■  (65) 

Then  substitute  Eqs.  (63),  (64),  and  (65)  into  Eq.  (55a).  On  using  the  first 
part  of  Eqs.  (3),  one  obtains  finally 


u  <x  ,  t) 

-v 


L 

x  =  ±l 


(66) 


where  w  (x  ,  t)  is  given  by 


w  (x,t)  f*  dp  e^H-Ve-vP^f 

A  "*■  0 


(p) 


Q  A*  M  f*  dP  e'lp(a-  5J  e-*,p\  ><p>  . 


(67) 


Clearly,  u  (x  ,  t)  of  Eqs.  (66)  and  (67)  is  a  sum  of  stratified  vorticity  modes 
u_  a(x, |  U  )  given  by  Eq.  (53a)  where  the  amplitudes  now  include  an  explicit  time 
factor.  It  is  seen  that  even  if  there  is  a  wind  which  drives  the  vorticity  modes, 
the  modes  with  smaller  pitch  or  wavelength  damp  out  much  more  readily  than  the 
longer  wavelengths.  It  was  our  hope,  when  we  first  investigated  this  solution, 
that  this  vortex  motion  would  have  some  of  the  characteristics  of  turbulence  in 
which  longer  wavelength  modes  feed  shorter  wavelength  modes.  However,  as  can 
be  seen,  stratified  vorticity  motion  does  not  have  this  property.  We  have  in¬ 
vestigated  other  exact  solutions  of  the  incompressible  Navier-Stokes  equation 
in  which  the  non-linear  term  v-V  v  plays  a  more  central  role.  Again,  however, 
these  special  solutions,  which  will  be  discussed  in  later  reports,  do  not  have  the 
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properties  of  turbulent  flow.  We  hope,  in  later  work,  to  discuss  more  general 
solutions  of  the  Navier -Stokes  equations  through  the  use  of  perturbation  theory. 

In  this  report,  one  of  the  solutions  of  the  linearized,  compressible  Navier- 
Stokes  equations,  which  we  shall  discuss  shortly,  has  some  of  the  properties  of 
turbulence.  However,  since  this  flow  is  longitudinal,  it  is  probably  not  what  is 
usually  thought  of  as  a  turbulent  flow. 

4.  2.  2.  '2  Stratified  Flows  in  Oceans  on  a  Rotating  Earth. 

We  will  now  consider  the  most  general  incompressible  viscous  flow  which 
can  occur  in  oceans  on  a  rotating  earth  subject  to  the  condition  that  the  velocity 
field  and  pressure  depend  on  the  vertical  coordinate  only.  It  will  be  seen  that 
sue1'  flows  are  stratified  vorticity  motions  of  the  type  just  discussed,  whose  time 
dependence  include  an  oscillatory  factor  which  depends  on  the  latitude  and  con¬ 
stitutes  a  sort  of  "tide". 

Following  Brunt  (1952),  the  equations  of  motion  of  a  fluid  on  the  rotating 
earth  expressed  in  terms  of  Cartesian  coordinates  in  which  the  x-axis  is  horizon¬ 
tal  and  points  east,  the  y  axis  is  horizontal  and  points  north,  and  the  z-axis  is 
vertical  is 

3v(x,t)  , 

—  +  I  v  (x ,  t)-  V  ]  v  (x ,  t)  +  —V  P(x  ,  t)  +  2  u  v  v  (x  ,  t) 

+  gk  -v  V2v  (x_,  t)  =  0  ,  (68) 


V  ■  v  <x  ,  t)  =0  , 

■v  "W 

where  «  is  the  angular  velocity  of  the  earth  about  the  North  Foie,  and  k  is  the  unit 
vector  in  the  direction  of  the  vertical.  Note  that  from  the  second  part  of  Eq.  (68), 
v  (x  ,  t)  can  be  expanded  in  terms  of  rotational  modes  as  follows: 

v(x,t)*V(t)  +  (2i)’3^  £  /dp  V  ■  (69) 

A.±l 


In  Eq.  (69),  Vt  i  is  a  thus  fax  undetermined  function  of  time  but  not  of  space. 
While  vectors  v^(x_,  t),  expanded  m  Eq.  (69),  satisfy  the  second  part  of  Eq.  (68), 
this  expansion  is  the  only  possible  expansion,  if  we  assume  that  we  are  working 
in  an  infinite  fluid,  that  is,  the  scale  of  variation  in  space  of  the  fluid  velocity  is 
small  compared  to  the  distance  cf  the  boundaries  of  the  fluid.  Whether  the  fluid 
is  of  finite  or  infinite  extent,  a  solution  wil  be  found  for  stratified  flow  which  is 
unique  in  the  case  that  xe  fluid  is  of  infinite  extent. 
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The  requirement  that  v  (^,  t)  depends  on  the  vertical  coordinate  z-x3  only 
leads  to  the  following  condition  on  gx(p  ,  t)  : 

gx(^,  t)  =6(px)6(py)gx(pz.t)  ,  (70) 

where  gs  (p  ,  t)  is  now  a  function  of  p  and  t  only. 

A  z  z  -3/2 

Then,  on  absorbing  a  factor  of  (2tt)  ’  into  the  function  gx  and  on  renaming 

the  variable  pz,  our  expression  for  v  (x_,  t)  becomes 


v  (x  ,  t)  =  V(t)  +  u  (z,  t)  , 


(7  1) 


where  u(z,t)  is  given  by 


/+0O 

dp  elp  gx(p,  t) 

•00 


(7  la) 


X=±l 


and 


Qx  =  Qx(k)  =  (2)‘1/2X(l,iX,  °)  .  (7 lb) 

It  is  clear  that  the  vector  u(z,  t)  is  the  sum  of  two  stratified  vorticity  modes 
with  the  direction  n  «  k. 

Note  that  uz  =  u3  =  0.  Furthermore,  the  non-linear  term  simplifies  in  an 
essential  way: 


3u  (z,  t) 

I V (5.' ‘>  •  1 1 t}  =  vx(t)  — 


(72) 


Let  6  be  the  co-latitude  measure  from  the  North  Pole.  Let  o>»|u>|be  the  a  - 
solute  value  of  the  angular  velocity  of  rotation  of  the  earth,  a«2u  cos  0,  0  =2u 

A 

sin  0,  V  (t)  be  the  vector  whose  components  are  [VX(D.  V  (t),  0),  and  W  (t)  be 
defined  by  W(t)  =  Vz(t). 

Finally,  the  dyadic  A  is  introduced  whose  components  are  given  by  the  matrix 
A  where  A  is  defined  below: 


/  0 


A  « 


-1  0 

0  0 


(73) 


0 


0 


0 


22 


The  x-  and  v-  components  of  Eq.  (68)  can  be  written  as 

A 


3V(t)  ~  3u(z,t)  3u  (z,  t) 

+  ffAV(t)  +  iHW(t)  +  - -  +  W(t)  - - 

+  aA  u  (z,  t)  -  y  (z,  t)  =  0  , 


(74) 


where  _i  is  the  unit  vector  in  the  x-direction:  _i  =  (1,  0,  0).  The  z-cornponent  of 
Eq.  (68)  is 


aw  (t) 
at 


0Vx<t)-0ux<z.t)+  j 


+  g  '  0  . 


(75) 


Equation  (74)  will  first  be  solved  by  substituting  Eq.  (71a)  for  u  (z,  t).  Note 
that  the  vi  ctors  gx  are  eigenvectors  of  the  dyadic  A  for  X  =  ±l: 

£3x  =  *iX2x  •  <76> 


Then  the  orthogonality  of  the  vectors  Q  x  is  used  and  a  Fourier  transformation 
is  taken  with  respect  to  z.  On  defining 

gx(o  “2x*  £(t)  •  (77) 

The  following  Equation  is  obtained  from  Eq.  (74): 

r  acx<t>  _  1/2  ] 

-gf - iA<*Gx(t)  +  (2)  ^fiW  (t)  J  6(p)  + 


3gx(p,  t)  „ 

— ^ - +  (-iAar  +  ipW(t) +ep  jgx(p,  t)  *  0  .  (78) 

From  the  theory  of  distributions,  the  coefficient  of  the  delta  “function  must 
vanish  as  must  the  second  line  on  the  right  of  Eq.  (78).  Thus  Gx(t)  is  given  in 
terms  of  a  differential  equation  as  dependingon  W(t).  Note  that 

'  £  Gx(t)gx. 

X»±  1 


V(t) 


(79) 
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On  integrating  the  differential  equation  for  G^(t)  and  then  on  substituting 
into  Eq.  (79),  one  obtains  the  following  relationships  between  the  components  of 
V(t>: 

V(t)  -  V{0)  cos  at  -  A  V (0)  sin  at  -  &  f  *  dt'  G(t-t')W(t')  , 

~  ~  ~~  Jq  ~- 

V1<t)  -  W(t)  ,  (80) 

where  G(t)  is  a  vector  Green's  function 

G(t)  *  (cos  at,  -  sin  at,  0)  .  (80s) 

A 

Thus  the  vector  V(t)  can  have  only  the  form  given  by  Eq.  (80),  where  V(0) 
and  W (t)  are  arbitrary. 

The  expression  u  (z,  t)  can  also  be  found  in  much  the  same  manner  as  in  the 
previous  section,  where  the  most  general  solution  for  the  stratified  vorticity 
mode  was  found.  The  differential  equation  for  gx(p,  t)  is  solved  by  the  usual 
elementary  methods,  and  the  reality  condition  is  used  on  u  (z,  t)  to  obtain  relations 
between  the  constants  of  integration.  Our  final  result  is 

u(z,t)«  £  wx(z-/*  dt' W(t'),  t).  (81) 

X.±i~A  Jo 


where 

wx(z,t)  ■  dP  eipZ  e  ^  ^(p) 

+  Q  *e'lX"t  dp  e"lpz  e"*'p2tfx*(p>  (81a) 

J  o 

Note  that  if  the  rotation  rate  of  the  earth  were  zero  or  if  we  were  at  the 
Equator,  the  velocity  u  (z,  t)  is  identical  to  the  stratified  vorticity  mode  Eqs.  (66) 
and  (67)  when  Let  us  now  take  f^(p),  which  in  general  is  arbitrary,  to  be 

given  by 

fx(p>  «  A6(p-k)  ,  (82) 

where  A  and  k  are  positive.  It  is  clear  that  the  general  w  x(z,  t)  is  a  super - 

_i  if 

position  of  such  particular  velocities.  On  absorbing  a  factor  (2)  '  into  A,  one 

obtains 


(83) 


w^x(z,t)  »  AX  co6fet+*kz)e  vk  \  w^fz.t)  *  -A  sin(at+Akz)e"1'11  X 
wXz<*'t>  ‘  0 

If  V(t)  is  picked  to  be  identical! _•  zero  and  one  of  the  chosen  to  be  equal 
to  zero,  then  the  velocity  is  given  by  Eq.  (83).  This  velocity  is  seen  to  be  a 
circularly-polarized  wave  travelling  in  the  negative  z-  direction  and  polarized  in 
the  direction  of  motion  if  X»  l,  and  travels  in  the  positive-  z  direction  and  is 
polarized  opposite  to  the  direction  of  propagation  when  A*-i.  If  the  viscosity  is 
not  zero,  this  wave  is  damped  in  time. 

It  is  now  easy  to  integrate  Eq.  (75)  for  the  pressure,  since  W (t),  Vx(t), 
ux(z,t)  are  known.  Let  us  consider  the  special  case  where  the  velocity  is  given 
by  Eq.  (83).  Then  the  pressure  is  given  by 

2 

P(z,  t)  -  PQ(t)  -  gz  +  £  Ap^Xe‘l'k  X  sin  ^  cos  (at+A  ^  .  (84) 

where  PQ(t)  is  the  pressure  at  z««0  and  is  an  arbitrary  function  of  the  time. 

It  is  perhaps  Interesting  to  note  that  Eqs-  (83)  and  (34)  provide  a  solution 
for  possible  motions  under  the  surface  of  a  "quiet"  ocean  in  which  the  surface  of 
the  ocean  is  taken  to  be  at  z«0.  Then  PQ(t)  is  the  atmospheric  pressure  at  the 
surface.  Such  solutions  are  the  simplest  vorticity  solutions  in  the  ocean  which  it 
la  possible  to  have. 

4.2  3  PROPAGATION  OF  SMALL  DISTURBANCES  IN  A  COMPRESSIBLE, 
VISCOUS  GAS 

This  section  treats  the  propagation  of  a  small  disturbance  in  a  compressible, 
viscous  gas.  It  will  be  assumed  that  the  fluid  flow  is  adiabatic  and  the  flow  is 
general  in  that  it  may  have  vorticity  as  well  as  irrotational  flow.  The  restric¬ 
tion  to  adiabatic  flow  is  not  essential,  for  equations  could  be  included  having 
energy  transport.  The  adiabatic  flow,  however,  is  simpler  and  illustrates  more 
directly  the  use  of  the  modes  which  have  been  introduced.  It  should  be  mentioned 
that  Morse  and  Ingard  (1968)  have  treated  closely  related  problems.  They  take 
into  account  heat  conduction.  They  also  Include  the  case  where  vorticity  is 
present,  but  they  separate  the  vorticity  by  using  well-known  vector  calculus 
relations  instead  of  using  the  curl  eigenfunctions,  as  in  this  report.  Further¬ 
more  they  are  Li  teres  ieu  primarily  in  time-harmonic  motions  and  in  solutions 
near  boundaries.  By  contrast,  our  interest  is  mainly  in  initial  value  problems 
solved  In  the  infinite  do  nain. 

Following  Landau  and  Llfshitz  (1959),  the  equations  of  motion  for  a  viscous 
fluid  are 


( V  (x  ,  t) •  V  J  v  (x,  t)  +£  V  P(x.t)  -^v2v(x,t) 

'■v  -V  fyo.  p  <V  P  <V 

t(W*)  £1£x.<*«‘>1  -  °* 

— fr —  +  (v(x  ,  t)-  V  Jp(x  ,  t)  *-p  (x,  t)V-  v  (x  ,  t)  ■  0  .  (85} 

C)\  «W  "W  <V  »V  ^ 

In  the  first  part  of  Eq.  (85),  p  is  the  usual  "dynamic"  viscosity  coefficient  and  t 
is  the  "second"  or  "bulk"  viscosity  coefficient. 

It  will  also  be  assumed  that  the  pressure  is  a  function  of  density  only: 

P  -  F(P)  •  (86) 


8v  (x ,  t> 

»v 


3p(x ,  t) 


Equation  (86)  holds  for  adiabatic  flow.  It  also  holds  for  isothermal  flow. 

We  will  now  investigate  how  a  small  disturbance  propagates.  This  disturbance 
is  regarded  as  toeing  a  disturbance  about  a  steady  state  whose  velocity,  density,  and 
pressure  are  constants  denoted  respectively  by  V,  p^,  and  Pq,  and  is  written  as: 

v  <x ,  t)  «  V  +  jj/x.  t) 

P(X.  t)  -  P0  +  P<X  ,  t) 

p(x,  t)  -  P0  +  p(x,t)  .  (87) 

In  Eq.  (87)  u .  P,  and  p  are  regarded  as  small  quantities  as  in  the  usual  perturbation 
approaches.  By  substituting  Eq.  (87)  in  Eq.  (85)  and  neglecting  products  of  the 

small  quantities,  the  linearized  equations  of  fluid  flow  are  obtained  below: 

» 

3u(x,t)  2  „  „ 

—"XT - +V?u(x,t)  +  V  p(x,  t)  -  vV  4u  (x  ,  t) 

Ol  "V/  ~  •'v  Pq  -W 

-4-(v+45)V[V-u(x,t)j  •  0  , 


&p(x,  t)  ^ 

—Jf~  +  £-£p< *■*>  +  P0l  I  '  0  • 

P(x  ,  t)  »  c2p(x  .  t)  . 

•>*> 


(88) 
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In  Eqi  (88),  we  use 


v  »  &i/ Pq)  *  kinematic  viscosity, 

f  -  <3C/4po),  c2  -  F'(p0)  .  (88a) 


In  Eq.  (88a),  the  prime  means  the  derivative  with  respect  to  the  argument. 
The  quantity  c  in  the  usual  treatments  of  linearized  compressible  flow  turns  out 
to  be  the  velocity  of  the  pound  waves.  It  will  play  the  same  role  in  the  present 
treatment.  The  first  two  parts  of  Eq.  (88)  are  regarded  as  being  a  pair  of 
simultaneous  linear  partial  differential  equations  for  the  unknown  quantities 

A  A 

u(x  ,  t)  and  p(x  ,  t).  The  last  part  of  Eq.  (88)  allows  one  to  obtain  P(x  ,  t)  from 

A.  A»  >V  V. 

P(x  .  t). 

From  Eq.  (18),  it  may  be  written  that 


u(x,  t) 


(89) 


where  the  summation  over  X  includes  all  three  values. 

It  is  clear  that  since  p  is  a  scalar,  it  should  be  expanded  as  an  ordinary 
Fourier  integral: 

p(x .  t)  -  (2  »f  3/2/ iji  elE. ‘  £.r{jj ,  t)  .  (90) 

Equations  (89)  and  (90)  will  be  used  in  the  first  two  parts  of  Eq.  (88),  and 
ordinary  differential  equations  will  be  obtained  in  the  variable  t  for  the  functions 
g^(p ,  t)  and  r(£,  t)  which  are  uncoupled  to  a  high  degree. 

On  substituting  Eqs.  (89)  and  (90)  into  the  second  part  of  Eq.  (88),  on  using 
Eq.  (12).  and  finally  on  taking  a  Fourier  transform  with  respect  to  x ,  there  is 
found  one  of  the  equations  we  wish  to  solve: 


®r(p ,  t) 

+  ’  iPpogo^E,,t>  *  0  •  <91> 

( In  Eq.  (91)  p  *  j  p  j ,  as  usual]  . 

Let  us  now  substitute  Eqs.  (89)  and  (90)  into  the  first  part  of  Eq.  (88),  and 
use  the  second  and  third  parts  of  Eq.  (12).  On  taking  the  Fourier  transformation 
with  respect  to  x,  one  obtains: 


j 


t 
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+  P 


(V  •  ^)gx(g_,t)  +i/p2gx(^,  t)j 
^  ir(p- l)  ‘  3"  <*/+4€>P«0<P  ■  t}j 


■  0  . 


(92) 


On  substituting  -pQQ(p  )  for  pin  Eq.  (92)  and  using  the  fact  that  the  vectors  Qx(p) 
are  linearly  independent,  one  obtains 


yg0(P*t) 

at 


+ 


4  2_ 


i(X.'S.)g0(S-'  +  3p  t2q(P  ■  tl  "  1  jT  Pr(P  *  l) 


■  0 


(93) 


where 

T  m„  +  §, 


and 


dgx(p  ,  t)  2 

— Til* —  +  i(V-p^)gx (p,  t)  +  rp  gx(p  ,  t)  *  0  (94) 

for  A  «  ±1  . 

Note  that  the  differential  equations  for  gx(p  ,  t)  uncouple  completely  for  the 
voi  tieity  modes  A  1.  Equations  (91)  and  (93)  couple  the  density  amplitude 
r(j)_,  t)  and  the  irrotational  velocity  amplitude  gQ  (p ,  t)  .  Thus  the  vurticity  dis¬ 
turbances  are  completely  uncoupled  from  the  pressure  and  density  disturbances. 

In  fact,  l-  q.  (94)  would  also  be  valid  fur  incompressible  flows,  if  perturbation 
theory  had  been  used  on  the  incompressible  Navier-Stnkes  equation.  In  a  sense, 
the  discussion  of  the  problem  of  small  disturbances  in  a  compressible  fluid  bears 
strong  resemblances  to  our  earlier  discussion  on  Maxwell's  Equations. 

Since  the  rotational  and  irrotational  components  of  u  (x  ,t)  uncouple,  it  is  con¬ 
venient  to  write 

u  (x^,  t)  «  u  (x,  t)  +  u,ir(x  ,  t)  ,  (95) 


where  u^  arid  m ^  are,  r<  spectivelv,  the  rotational  and  irrotational  parts  of  u 
and  are  thus  given  hy 


(95a) 


u  ir(x.t)  -  /jd  x  ofc|E.)B0<E.*t>  ■ 

In  Kq.  (05),  the  real  part  of  the  expressions  are  taken  because,  as  explained 
in  connection  with  1'q.  (13c),  the  reality  conditions  on  the  victors  lead  to  con¬ 
ditions  on  the  amplitudes  g^(p,  t)  which  allow  us  to  assume  that  the  imaginary  part 
of  the  expansions  are  zero. 

liquation  (94)  is  easily  solved.  The  solution  is,  in  fact, 

*  6a<P^*P  l-‘(V-p)t-  rp2t),  (99) 


where  g^(p)  ■  gA(p  .  0)  is  the  constant  of  integration. 

We  have  immediately 

u  t.(x  ,  t)  »  X!  w^(x-Vt,  t),  (97) 


where 


w  A(x  ,  t) 


*  A(i|P)eA<P>,'rp  1 


(97  a) 


it  is  clear  that  u^  is  a  superposition  of  vorticity  modes  which  are  being  car¬ 
ried  along  with  the  steady  wind  who^c  velocity  is  V  and  whose  amplitude  is  being 
attenuated  with  the  time  factor  e'l’P  l.  Hence  the  vorticity  modes  with  shorter 
wavelength  damp  out  at  a  much  faster  rate  than  those  with  a  longer  wavelength. 
That  is,  the  small  scale  vorticity  modes  are  damped  out  faster  than  the  larger 
scale  modes.  Thus,  the  rotational  velocity  uf  does  not  at  all  behave  like  tur¬ 
bulence  in  this  respect. 

We  shall  now  show  that  the  irrotational  velocity  u  .j,  and  the  corresponding 
density  increment  p  can  be  written  as  the  sum  of  two  parts:  one  which  can  be 
interpreted  as  a  damped  sound  wave,  and  the  other  which  is  a  damped  standing 
wave.  The  damping  of  the  sound  wave  is  such  that  the  longer  wavelengths  persist 
longer,  while  for  the  standing  wave  the  shorter  wavelen  ths  maintain  themselves 
for  a  longer  time.  Surprisingly,  then,  this  standing  irrotational  wave  has  some 
of  the  characteristics  of  turbulent  flow  in  that  small  scale  phenomena  persist 
over  large  scale  ones. 
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Differential  E’qs.  (91)  and  (93)  will  now  be  solved.  Let  R(^,t)  and  G(jo,t)  be 
defined  by 


r(p,  t)  *  e  *E.  X-^TKjp,  t).  gg(£.  t)  «  e  ^  ^.^(p,  t)  .  (98) 


The  differential  Eqs.  (91)  and  (93)  are  now  differential  equations  for  R  and  G: 


3fi(p.t) 

— 3? - ipopG(p  ■ l)  “ 0  - 


ac(p,t)  4  2 

- - +3.Tp*G(p,t)  -  i£-pH<p,t)  «  0  .  (99) 

dt  J  pO 

The  usual  techniques  are  used  to  solve  these  coupled  differential  equations. 
The  second  part  of  Eq.  (99)  is  differentiated  with  respect  to  t,  the  0H(j^,  t)/i)t  is 
eliminated  through  the  use  of  the  first  part  of  Eq.  (98),  and  the  following  second 
order  differentia!  equation  is  thereby  obtained  for  G: 

A(|  3G(p.t) 

-  -f  •  .c  Tp“  — q?2 —  +  c  p  G(p,t)  =  0,  (Kill) 

«—  »*  "w 


The  general  solution  of  this  equation  is  of  the  form 

G(£_.  t)  *  Cl+(p)exp(v/+tJ  +  G_(p)cxp(w_tJ  ,  (101) 

where  w  an-  tie  two  -  outs  of  the  quadratic  equation 

2  4  >  ••  v 

w  +  3  TP''-  J  ‘  0  ,  (101a) 

and  G±(p)  arc  i  ei.su.nts  of  integration.  The  plus  and  minus  signs  refer  to  the 
signs  in  front  of  the  disci  in. intint  of  the  quadratic  equation.  Of  course,  w 
depends  on  p. 

H<  cause  of  tin  radical!',  dilfercnt  chai'acter  of  w  4  for  different  regions  in 
p -space,  it  is  coma  incut  to  ■  ivide  p -space  into  two  domains.  A  critical  length 
L^  is  introduced  ami  a  rori'csp. Hiding  <'ritical  wave  muoln  i  U  by 

k  ”  f~  1  <3c/2t)  .  (102) 


Call  P  the  sphere  in  p  -space  for  which  p  <  k,  where  the  letter  P  stands  for 
"Propagating".  The  remainder  of  j^-space,  that  is  all  of  p  -space  other  than  the 
sphere  P,  will  be  denoted  by  N,  which  stands  for  "Non-propagating".  Now  the 
results  can  be  combined  to  give  gg(p,t)  and  r(p  ,  t)  in  a  physically  significant 
marner. 


For  p  in  P  : 

gfi(p  ,  t)  *  exp'  -ip  •  Vtjexp|  -cp(p/k.'t]  |c+(p  >exp|  icp(l-(p/k)2)1/|,2t] 

+  G_(p  )exp[  -icpO-fp/k)2)1/2^  |  , 

Pn  ( 

r(p  ,  t)  =  -i-j-  erp[  -ip  •  Vt]  exp[  -cp(p/k)t]  |  G+(p  ) 

v  [  (p/k)+in-(p/k)2)1/l2]exp[icp(l-<p/k)2)1//2tj 

+  G_(p  )[  (p/kJ-id-Cp/k)2)1^2]  exp[  -icp(l- (p/k)2)1/,2t  J  j-  .  (103) 

For  p  in  N: 

g0(EL't)  =  exPt-iE.-VJ  {G+(p)expI-CpF+(p/k)tj 

+  G_(p  )exp[  -cpF_(p/k)tJ  }  , 

~  ' 

Pn 

r(p.t)  =  -i~  exp[  -ip  -  V  J  |  G+(p  )F_(p/k)exp[  ~cpF+(p/ktJ 

+  G_(p  )F+(p/k)exp[ -cpF_(p/k)t]  |  ,  (104) 

where  the  function  F_^(x)  which  appears  in  Fq.  (104)  is  defined  by 

F+(x)  =  x  -  (x2-  1) 1^2,  F_(x>  =  x  +  <x2-  l)^^2  .  (104a) 

On  using  Eqs.  (103)  and  (104)  in  Fqs.  (95a)  and  (90)  it  is  seen  that  u;j  and 
p  each  can  be  broken  up  into  two  parts  corresponding  to  the  spaces  P  and  N.  It 
is  also  seen  that  the  parts  of  u  ^  and  p  that  correspond  to  P  space  are  damped 
moving  waves  which,  >r  sufficiently  small  p,  are  undamped  and  travel  with  the 
velocity  c  with  respect  to  the  constant  wind  which  has  the  velocity  V.  The  shorter 
wavelengths  damp  out  faster  than  the  longer  wavelengths  (with  a  different  damping 
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factor  from  that  for  rotational  flow),  so  that  the  longer  wavelengths  persist  for  a 
longer  time.  Thus  the  contribution  of  the  space  P  may  be  identified  as  being 
sound  waves. 

The  contribution  to  u.  and  p  from  the  space  N  does  not  at  all  behave  like 

'■wlx 

sound  waves.  These  motions  are,  in  a  sense,  "overdamped". 

Because  of  the  nature  of  the  function  F  which  appears  in  Eq.  (104),  it  is 
seen  that  the  shorter  wavelength  contributions  damp  out  more  slowly  than  the 
longer  wavelength  contributions.  Thus  small  scale  phenomena  persist  longer 
than  the  large  scale  phenomena.  In  turbulent  flow,  a  similar  pattern  is  observed, 
though  usually  it  is  assumed  that  turbulent  flow  is  rotational,  rather  than  ir- 
rotational.  Our  results  suggest  that  perhaps  some  of  the  phenomena  of  turbulent 
flow  are  related  to  the  irrotafional  components  of  the  How  corresponding  to  the 
space  N.  Perhaps,  also.  Inclusion  of  the  non-linear  terms  in  the  equations  of 
flow  would  induce  a  similar  vehaviour  in  the  rotational  flow.  However,  such 
speculation  is  beyond  the  intent  of  this  report. 

Without  going  into  detail,  note  that  all  the  constants  of  integration  g^(p  )  and 
G±(p  )  are  ur.iquely  determined  when  u  and  p  are  given  at  time  t=0.  Hence  the 
initial  value  problem  has  been  solved  in  all  generality  for  the  propagation  of  small 
disturbances  in  a  compressible  viscous  fluid. 


3.  H  HTIIHR  PROmtTIKS  Oh  THE  \  KETOH5  tf^)  AM) 
(iENEIUEIZEl)  sI  KKUiE  I! XKMOMES 


Having  given  some  applications  of  the  eigenfunctions  of  the  curl  operator,  we 
shall  go  into  further  properties  of  the  vectors  Q  >(o),  primarily  to  obtain  the  ro- 
tation  property  of  Eq.  (9)  and  to  restate  the  Helmholtz  theorem  in  terms  of 
spherical  coordinates.  Up  to  now,  it  will  be  noted  that  the  cartesian  coordinates 
are  emphasized  over  other  coordinate  systems. 

It  will  be  necessary  to  go  into  some  properties  of  the  representation  of  the 
rotation  group.  Indeed,  as  will  be  shown  in  the  Appendix  A,  rotation  group  con¬ 
siderations  have  motivated  our  introduction  of  the  Helmholtz  theorem. 

The  irreducible  representations  of  the  generators  of  the  rotation  group  are 
three  matrices  (i*l,  2,  3)  which  satisfy  the  commutation  rules 


[s!(j)*  s2(j)]  *iS30)  <°yc-  > 


(105) 


where  j  characterizes  the  irreducible  representation  and  is  a  positive  integer  or 
half-odd  integer.  The  matrices  have  2j+l  by  2j+l  in  dimension.  The  elements 
are  denoted  by  (j,  m| Sj j,  m')  where  m,  m 1  take  on  the  values  -j,  -j+i,  . . .  ,  J  - 1,  j. 
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The  usual  form  for  the  elements  are  (for  example,  see  Edmonds,  1957). 
(j.m|s3|j,m->  =  mAmm,  , 

0,m|Sj|j,  m')  ±  i(j,m|S2|j,  m')  =  (  (jTm'Xjim'+l)]  rai±1’  (106) 

The  column  vectors  upon  which  the  matrices  act  will  have  the  components 
labelled  by  m,  where  the  top  component  will  be  labelled  by  m=j,  the  second 
component  by  m=j-l  and  so  on  until  the  bottom  component  which  is  labelled  by 
m  =  -j.  Correspondingly,  the  first  row  of  the  matrices  S.^  will  be  labelled  by 
m=j,  the  second  by  m=j-l,  etc.  A  similar  labelling  convention  is  used  for  the 
columns  of  these  matrices. 

Let  0  =  (0^,  6^,  0g)  be  any  real  vector.  Then  exp  [  i6^  S  is  an  irreducible 
representation  of  the  rotation  matrix  R(0)  of  Eqs.  (7)  and  (8).  (We  use  0-S^  = 
S6,S.^).  Moses  (1966b)  has  given  the  matrices  exp[i0-S^']  for  all  0,  It  will  be 

[  1  i  /  o  \  ~  ^ 

necessary  to  introduce  the  Jacobi  polynomial  (*)  •  Its  properties  are  given 

by  Szego  (1959),  for  example.  We  define  the  function  S(j,  m,  m ’,  x)  by 


SQ.m.m'.x)  =  pto-m'.m+m*)  (x)  . 

From  Rodrigues'  formula  for  the  Jacobi  polynomials  (Szego,  1959) 
S(j,m,m',x).(-l)j-m  (i-x)”(rn”rn>(l+x)'(rn+rn ') 

[a^-  d+x)j+m'l. 

dxJ  L  J 

Moses  (1965)  has  given  an  alternative  expression  for  S(j,  m,  m  x): 


S(j,  m,  m',  x) 


(-l)j+m' 


2™-j 

(j+m)! 


d^m 

dxj+m 


(107) 


(1C8) 


(109) 


In  fact,  one  may  use  either  Eq.  (108)  or  (109)  to  define  S(j,  m,  m ',  x).  Let  us 
denote  the  matrix  elements  of  exp  (iQ’S^J  by  (j,  m |exp[  i0  •  S  j|  j,  m!).  Then 

(j,  m|exp|  ii9  •  Sj|  j,  m ')  »  ( (j-m)l  /(j-m1) !  ]  1/,2[  (j+m)!  /(j+m')!  ] 


.  9sm-m'  /  ei+ie2\m*raV  8X.93  .  ,  .  (110) 

c(sm  g-)  (  — g - J  (  cos  y  +  l7f  siny  )  S{j,m,m',z)  , 
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where 


,  r 

/Ml  /M 

1*  2  =  [ 

1  *(t)  j  cos0+[e)  • 

(110a) 

We  have  previously  introduced  "generalized  surface  harmonics"  (Moses, 
1967 ).  The  generalized  surface  harmonics  Yj11'111  are  defined  for  0<6<n 
and  O<0<2ir  by 

Ym'm'  (f),4>)  =  (-l)m'm'(l/2)m  +  1[(2j+l)/ir]  1/2[(j-m)!/(j-mt)!  j  1/2 

xKi+m)!/(,i+m',!]l/2ei<m-m'^[sin  e]™-"1’ 


x  (  '+cos  0Jm  S(j,m,m',cos  fl].  (Ill) 

These  functions  have  the  following  properties: 

Y™'0<6.d)  *  Yjm(0,tf)  ,  (112) 

where  Yjm(6.  d)  are  the  usual  surface  harmonics  in  the  notation  of  Edmonds. 

f  2  V  f  TT  , 

Jfo  dd  d0  sin  0  Ypn  <0,d)Y™  'n(0,  d)  *  4.  jj,i  »  (113) 

00  J 

E  £  Y"ljn(6,  d)Y™,n  (O'.d')sin  O'  -  6(O-0')6(d-d')  ,  (114) 

j=|n|  m-j 

j 

£  Y™,n  (O.d)Y?,’n'*(0.d>  =  !(2j+l)/4ir)6n  n,  ,  (115) 

m=-j 

j 

£  Y™'n(O,d)Y.rnl'nt(0,d)  *  [(2j+l)/4Trj6m<m 
n=-j 


(116) 


4 
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Y™'n  (6,4)  =  (-l)n‘mY  n,m*(0,$>  ,  (117) 

J  J 

j+k 

Y^'n'  (0,  (6,4,)  *  £  |  ((2j  +  l)(2k+l))/(4-n-(2J  +  l)j|  1/2 


x  (k,  n',  j,  m'|  k,  j,  J,  n'+m '  )(k,  n,  j,  m  |k,  j,  J,  n+m ) 
v  Y™+n'rn'+n’  (0,4,).  (118) 

In  Eq.  (118),  we  have  used  the  Clebsch-Gordan  coefficients  in  the  notation  used 
by  Edmonds. 

The  matrix  elements  of  exp[i£-S  for  the  case  6^*  0  is  conveniently  ex¬ 
pressed  in  terms  of  the  generalized  surface  harmonics.  Using  polar  coordinates 
to  describe  0  which  has  only  x  and  y  components,  0  and  4  are  defined  by 

6  *  1 0  j  6  =  0(cos  4>,  sin  <t>,  0)  .  (13.9) 

Then 

(j,m|exp[i0-S  jjj.m’)  *  (-i)m~m'  (  4ir/(2j+l>]  i/2Ynr.,  m  {Q  ^  _  (120) 

One  of  the  more  interesting  properties  of  the  generalized  surface  harmonics 
is  their  behaviour  under  rotations.  Let  the  angles  0,  <t>  determine  a  unit  vector 

Z  by 

7)  »  (sin  0  cob  4,  sin  0  sin  4>,  cos  0)  .  (121) 

Let  17  '  be  the  vector  obtained  from  r)  by  means  of  the  rotation  R(-fi  ).  That 
is 

n '  -  R(-«)n  .  (122) 

where  the  elements  of  the  rotation  matrix  arc  given  by  E.  (8).  Furthermore, 
let  O'  and  4>'  be  the  angles  which  give  rj'  through  Eq.  (121).  Then 


i 


I 


i 

i 
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j 

Y™n  (0\  ^')  »  exp  [  -2in^(£2  , )]  £  y"1'*  n(0,  *)(j,  m'|exp[ifi-  S]|j,  m)  ,  (123) 

m’«-j 


where  <t>(n,r))  is  given  by  Eq.  (9b)  . 


It  will  now  be  convenient  to  express  the  matrix  R(0)  whose  elements  are  given 
by  Eq.  (8  )  as 


R(£)  =  exp[i£  •  K]  , 


(124) 


3 

where  £•  K  *  £  0.K.,  where  K.  are  three  matrices  given  by 


It  is  a  well-known  theorem  that  Eq.  (124)  is  valid,  and  can  be  verified  by 

3  2 

direct  expansion  of  the  exponential  using  the  fact  that  ( 0 ■  K)  “0  (0  •  K)  . 

Now  the  matrices  Kj  are  unitarily  equivalent  to  the  matrices  S.^: 

K.  =  VSi(1)V+  (126) 


where  the  dagger  means  hermitian  adjoint  and  the  matrix  V  is  given  by 


V  =  (2)‘i 


/2 


It  is  clear  that 


(127) 


V1*  vt 


(128) 


u. 


ii  1!; 

[$•  I 


b 
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Let  us  construct  the  matrix  r  as  follows  for  any  unit  vector  r?  : 

r-  V  exp[iu-S<1)]  =  exp[iu.K|  V  ,  (129) 

where  u  is  related  to  rj  by  means  of 

d(sin  <{>,  -cos  <t>,  0)  ,  (130) 

where  >)  is  given  in  terms  of  9  and  $  by  Eq.  (121). 

Then  the  vector  Q  j(rj )  has  its  components  formed  from  the  components  of  the 
first  column  in  I’,  ihe  top  component  of  the  column  vector  being  the  x -component 
of  Q ,(i>),  the  second  component  being  the  y-  component,  and  the  third  component 
being  the  i -component.  QqOi  )  and  Q.  j(l )  are  formed  similarly  from  the  second 
and  third  columns  of  V  . 

The  rotation  properties  of  Q^(rj)  given  by  Eq.  (9)  then  follow  from  Eqs.  (120) 
and  (123).  Let  V  ,,  V_,  V  ,  be  vectors  constructed  from  the  first,  second,  and 

^•1  '•v  U  -w  ”  1 

third  columns  of  the  matrix  V  of  Eq.  (127).  That  is 
vx  ■  Cl)'1/2  <x,i,  o)  .  for  X  »±l  , 

V  Q  =  (0,  0,-1)  .  (131) 

Then 


Q  K(r] )  «  (4»/3)l/2£v  .  (132) 

6.  THE  HELMHOLTZ  THEOREM  IN  TERMS  OE  SPHERIC.  A  I,  COORDINATES 

This  section  gives  the  general  vector  u(x)  in  terms  of  spherical  coordinates 
in  such  a  way  that  the  Helmholtz  decomposition  u(x)  =  S  u_^(x  )  is  preserved. 

The  components  of  u^(x)  will  also  be  given  in  terms  of  spherical  coordinates. 
Toward  this  end,  it  will  be  useful  to  introduce  the  vector  spherical  harmonics. 

6.1  Vector  Spherical  Harmonics 

The  vector  spherical  harmonics  are  defined  in  terms  of  our  notation  (Blatt 
and  Weisskopf ,  1952)  by 


ijLM 


(0. 1)) 


£  Vx  Y1"'°(0,4>ML.m.  1,*  |  L,  l.J.M) 
m,  X 


(133) 


4 


t 


>-*•»!  41  IsajrthHK’B  MMlUVNIlit".-  -Ii 
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where  is  given  by  Eq.  (131)  and  the  summation  is  over  all  X.  Instead  of 
giving  the  properties  of  the  vector  spherical  harmonics,  we  refer  the  reader  to 
Blatt  and  Weisskopf.  We  shall  be  interested  in  the  components  of  Y  TT  <g,  ij>)  in 
spherical  coordinates.  Let^,  j  ,  k  be  the  three  unit  vectors  in  the  x,  y,  z  direc- 
tions,  respectively.  Let  a_^  be  the  unit  vectors  along  the  r-direction, 

O-direction,  and  $ -direction  respectively.  Then,  as  is  well  known. 


a_r  *  sin  0  cos  i  i  +  6in  6  sin  <t>  j  +  cos  0k, 
a  q  *  cos  0  cos  <j>  ^  +  cos  8  sin  <t>  j  -  sin  0  k  , 
*  -sin  <t>  i  +  cos  <t>  j  . 

The  following  can  be  shown: 

ir‘XjJM<0'^“  °* 

ir  •  It,  j-  1,  M<9'  ♦>  -  ( J1/2y}*'  °(0,  *)  , 


~r"  ~  J,  J+l,  M 


(0,  4)  *  -(  <J+1)/(2J+1)  ]X/2y“*  °(0.  «) 


sin  9  Id' 


M 


|J(J+l)j 


-JTi 


sin  9S_0-Ij,  J-1.M(0*^)  "  2l+TI(J'M+1){J+M+1)/J(2J+3>)  1/,Z 


x  Y  .  ' ,  (0,  4)  - 


■  I  <J-M)(J+M)/J(2J-1)]  »/aY“.>.*>  , 


J+l  2J+T 

sin  9i0'lj,  J+i,M(9,0>  *  27tTi(J-M+1>(J+M+1)/(J+l)(2J+3,]1/Z 
x  yJ^O,*)-  |j^[(J-M)(J+M)/(J+l)(2J-l)l  1/2y5**1°  (©»  . 

sin  6Z,<t,‘Z "  'l{uTtJ-M+l)(J+M+l)/(J+l)(2J+l)(2J+3)J1/2 
x  Y^[°(0,4)  -  !<J+1)(J-M){J+M)/J(2J-1K2J+1)]1/2Y^_'°(0,  «)|  , 


+  JWHH  ;W-'  :■ 


(134) 


(135) 


(136) 


/ 
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fc 


ain  0  a . 


M 


fJ(2J+l» 


T7?Yf  V'*> 


M 


[  (J+1)(2J+1)| 


1/2 


(137) 


In  Eqa.  (135)  through  (137),  one  must  recollect  that  Y™’  0  (0,  <t> )  are  just  the 
usual  surface  harmonics.  See  Eq.  (112).  The  following  sketches  the  derivation 
of  Eqs.  (135)  through  (137).  Note  that: 

a  •  V  „  ■  -  cos  6  , 

a  -Vv  -  (2)'1/2  sin  9  el  *  ,  for  1*11.  (138) 

~s,r 


"  Bin  6  * 

a^-V^  «  (2)"^2A  cos  0e'^  ,  for  A  ■  ±1  .  (139) 

a.-v.  ■  0, 

a^-V^  *  i(2)*1^2  e^  ,  for  A»±l  .  (140) 

Also 

cos  9  «  (4Tr/3)1/,2Y^°(0,*)  f 
sin  9  -  -(8Tr/3)^2Yj’O(0.<|))  , 

sin  9  e'1*  «  (8t/3)1/2Y1'1*  °(8.«  .  (141) 

After  multiplying  through  by  the  appropriate  unit  vectors  on  both  sides  of 
Eq.  (133)  and  also  multiplying  through  by  sin  9  to  obtain  Eqs.  (136)  and  (137),  sub¬ 
stitute  the  appropriate  expressions  from  Eqs.  (138)  through  (141)  on  the  right  hand 
side,  use  Equation  (118)  to  reduce  the  product  of  surface  harmonics,  and  finally 
use  the  explicit  forms  for  the  vector  coupling  coefficients  (L,  m,  1,A|  l,  1,  J,  M). 


6.2  Expansion  of  Vector  Fields  in  Teras  of  the  Irreducible 
Representations  of  the  Rotation  Group 

Let  us  now  return  to  the  problem  of  expanding  an  arbitrary  vector  u(x)  in 
terms  of  spherical  coordinates.  The  starting  point  is  the  expansion  of  Eqs.  (13) 
and  (13a).  Define  the  function  g(p,j,m,A)  by 


(142) 


g(p,  j,m,X) 


sin  8  dfl  Y1 


m,  X* 


<®.  ^>gx  <£>  • 


In  Eq.  (142),  8,$  is  used  as  the  polar  angles  of  jp;  it  is  desired  to  reserve 
8,  $  to  denote  the  polar  angles  of  _x.  Strictly  speaking,  g(p,  j.m.X  )  is  defined  for 
j  >  |x|  only,  but  g(p,  j,  m,  >-)  will  be  defined  for  j  »  -1,  0,  1,  2, . . .  with  the  under¬ 
standing  that 

g(p,  0,  0,  ±1)  ■  0  ,  g(p,  -l,m,A)  *  0  (142a) 


for  the  sake  of  minimizing  the  difference  in  appearance  of  future  expressions 
between  the  cases  for  which  X  «  0  and  X  ■  ±  1. 

From  the  completeness  relation  for  the  generalized  surface  harmonics, 

Eq.  (114),  gx<£_)  can  be  found  from  g(p,  j.m.X)  through 

°°  j 

*  X!  £  Yjn’  •  <J43> 

j»0  m»-j 


Thus  one  rnay  express  the  vectors  u^(x  ^  aB 

V1/2 

J.m  P 

,m,X 


«x(x  )  =  <60' ^  Y^x‘(e,*)v^ 


v  Yp'<0.*)g(p,  j.m.X)  , 


(144) 


where  Eq.  (132)  has  been  used  ior  Q  \ln  )  • 

It  is  clear  that  from  the  construction  and  from  Eq.  (12)  that  v  xu.lx  )  has 

•V  -VA  'Nr 

the  same  expansion  as  Eq.  (144)  except  that  g(p,j,  m,X)  is  replaced  by 
Xpg(p,  j<  m,  X).  Furthermore,  from  Eqs.  (12)  and  (13a),  it  will  be  shown  that 


»  k 

Z ‘Lo{*>  *  -(2/,r)1/2i  £  £  Y™'°(e,^)(i)K 

k“0  in'-k 


dP  P  Jk(pr)g(p,  k.m,  0) 


(145) 
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where  r  *  !  x  |  and  0,  < p  are  the  polar  angles  of  x  .  In  Eq.  (145),  the  functions 
Jk(x)  are  the  usual  spherical  Bessel  functions.  To  derive  Eq.  (145),  Eq.  (12) 
is  used  in  Eq.  (13a)  for  A«o.  But  a  well-known  expansion  in  terms  of  our 
notation  for  elE-’~  in 

oo  q 

elE.'I  »  4*  £  £  <i)qjq(priY^'°(0.*)Y^*0i,(e,J)  .  (146) 

q»0  r«-q 

Equations  (145)  is  obtained  by  substituting  this  expression  for  e'L  then 

A  A 

integrating  over  0,4  and  finally  using  the  orthogonality  relation  for  the  surface 
harmonics  of  Eq.  (113). 

rhe  functions  g(p,  j,  m,A  )  have  interesting  properties  under  rotations  of 
coordinates.  Under  the  rotation  of  axis  described  by  the  vector  0  ( 0  is  per¬ 
fectly  general  and  has  no  connection  with  the  polar  angles  0,  4),  the  vector 
u^(x)  goes  over  into  uA'(x)  where  u^'(x)  is  given  by  Eq.  (17).  Similarly,  the 
ami-litudes  g^(g^)  go  over  into  g^'(]^),  where  gx'fgj  is  given  by  Eq.  (18).  It  has 
been  shown  by  Moses  (1967,  1970)  that  the  relations  Eqs.  (142)  and  (143)  which 
relate  g(p,j,m.  A)  and  g^<£)  lead  to  the  following  relation  between  the  transformed 
coefficients  g'(p.j,rn.  A)  and  the  coefficients  in  the  original  frame  g(p,  j,m.  A): 

J 

g’(p,j.  m.  A)  a  V  (j,  m  exp  ( i0  ■  S  j  j,  rn')g(p,  j,  in'.  A)  ,  (147) 

m 1  *-j 

where  (j,  in |exp[  iO  ■  S  j|j,m')  is  given  by  Eq.  (110).  Thus  the  coefficients 
g(p.J.m,A)  transform  under  the  irreducible  representations  of  the  rotation  group 
without  any  mixing  of  the  variables  A.  This  expansion  confirms  the  rotational 
invariance  of  the  Helmholtz  decomposition  of  this  report. 

Equation  (144)  will  now  be  simplified  through  the  use  of  vector  spherical 
harmonics.  The  right  hand  side  of  Eq.  (146)  is  substituted  for  el6,’ ,  the 
product  of  Y^*  *■*  (0,  4)  and  Yj*  (0,4)  is  reduced  through  the  use  of  Eq.  (118), 
the  orthogonality  relation  Eq.  (113)  is  used  on  integrating  over  the  variables 

A  A 

0, 4  of  p,  the  definition  of  the  vector  spherical  harmonics  Eq.  (133)  is  used,  and 
finally  the  known  values  for  the  remaining  vector  coupling  coefficients 
(L,  0,  1.  A|  L,  1,  J,  A)  are  substituted  to  obtain  the  following: 
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oo  k 

,  \  ,wv-l/2  v-  v"  r\*  -XY,  .  (9*4)  A.  x  , 

ux(x)  *  -<ff)  '  2-  2^  (i)  ^k.  k,m  k,  m/,| 


(r) 


k“l  m*-k 


+  Ik/(2k+1)]i/2Yk>k+lam(6.«)Ak#mX  k+1(r) 

-i|  (k+l)/(2k+l)J  1/2Yk>  k.u  m(0.  *,Ak<  (  k.  ,(r)  j  , 

for  X  ■  ±  X,  and 

«o  k 

u  Q(x)  =  (2/*)1/2  £  j;  (Uk+1  j[(k+l)/(2k+l)]1/2Ykjk+1>m(e.1j) 

k*0  m»-k 

v  Ak,m,0.k+l<r>  ^k/(2k+l)il/2Yk 'k.litn(0.4) 


v  Ak,m,0,k-l(r) 


(148) 


where 


/WT: 

P2dp  jk(pr)g(p.j..n,X)  . 


(149) 


The  vector  potential  Ax(x_)  of  Eq.  (15)  for  X  *  ±  1  has  the  same  expansion 
as  the  first  of  L'q.  (148)  with  g(p,j,m,X)  replaced  by  (X/p)g(p,  j,  m,  X).  The 
scalar  potential  V(x_)  of  Eq.  (15)  is  given  by 


V<x)  *  <2/r)l/2  53  £  Y™’0<«,C>)  J  pdpjh(pr)(i)k+1g(p,l<.rn,0).  (150) 


k«0  m»-k 


Though  the  expansions  of  Kq.  (14  8)  are  very  useful,  for  many  applications 
one  wishes  the  components  of  ux(x)  in  spherical  coordinates.  Hence  ar-ux(x)  , 
sin  0  if"  X(x_)  and  sin  8  w ill  be  given  in  terms  of  g(p,  j,m,X).  One 

uses  Eqs.  fl35)  through  (137)  and  sees  that  the  components  of  u^bO  are  sums 
oi  series  in  Y™‘  .  In  order  that  the  components  be  a  series  in  Y™’  ® 

only  (that  is,  a  series  in  surface  harmonics),  k±l  is  replaced  by  k  in  sum¬ 
mations  and  the  series  is  then  expressed  in  terms  of  the  spherical  Bessel 
functions  Jk>  jk±T^k±2  *  functions  can  ^  expressed  in  terms 

of  Jk  and  derivatives  of  Jk  by  using  the  well-known  relations,  recursively,  if 
necessary: 
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few* 


=jk-i(x) 


7jk(x)  -  ^jk(x)  ■  jv+,(x>  . 


Jk+1 


(151) 


The  differential  equation  for  the  spherical  Bessel  functions  is  also  used,  namely 


~2  jk(x>  +  2x  ~  j.(x)  +  {xZ-k(k+l)]  j,(x)  =  ° 
dx  '  R 


(152) 


>Y  '  .  : 

M  "•  , 

\J  '>■ 
v;v  ■ 


<  mm 

K  Xi 

I ’  jit 

*j|P 

jr 

*,4$| 


m 

m 

: 

■*  > 

'?  ;?  « 

>  * 


,,4 

v  *:*"■ 
I:-  jL 

ft- ! 

i 


Though  the  details  of  the  calculation  are  tedious,  one  can  give  the  final  ansv  -rs 
in  a  straightforward  manner. 


oo  k 


lr  ■  Wl/2  jr  £  £  U)k+1  lk(k+l)J  1/2Y™'  °(«. 


#> 


k»0  m“-k 


fQ  p  dp  ik(Pr)8<P-  k»  m.x  )  , 


(153) 


fo.  a  *  ±  1,  and 


oo  k 


a.r-lic(x)  •  (2/^) 1/2  £  £  (i)k+1  Y“'°(9.*) 

k“0  m*-k 


J“  2 

P  jk(pr)g(p,k,  m,0) 


(154) 


In  E<;.  051 the  prime  means  derivative  w»,u  respect  to  the  argument  of  the 
function  Similarly,  later,  the  double  prime  means  that  a  second  derivative  has 
beer  taken  with  respect  to  the  argument. 


oo  k 

o-.re  ^-UX(X  )«-(»)' 1/2  £  £  (i)k  v«n,0  (  ^ 

k“0  m«-k 


+  ((k-i)fk-m)ik-Hr!,Vk?2k-l><2k+l)J1^2l  1-  Jg(p,  k-l,m,X) 


•  T» 
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-  (fk+2)(k-m  +  l)(k+m+l)/{k+l)(2k+l>(2k+3)j1/2^  -  k(k+l)  ^ 

x  g<p,k+l,m,X)J 

~  T  fQ  p  dP  £l(k-l)k(k-m)(k+m)/(2k-l)(2k+l)j  */2 

Y  C(p i  k-1,  m.X) 

*  I  (k+l)(k+2)(k-m  +  l)(k+rri+l)/(2k+l)(2k+3)|  */2 

l,m,X)jj, 


x  g(0,  k+1. 


(155) 


for  X  =  ±  i  . 


In  Eq.  (155,  and  later.  g(p,  0.  0.  ij)  -  g(p,-l,m,A)  =  0.  An  alternative  ex¬ 
pression  for  sir  0  a^  •  uA(x  )  for  X  =  ±1  is  given  by 


«J  k 


sin0lo-iL\<V  *  -  E  E  (i)kY^'V. 

k=0  m*-k 


mX 


f0  p  dP  Jk(pr)  [^k'1)(k*n>)(k+m)/k(2k-l)<2k+l)J  X/2 
g(p,  k- 1,  m,  X)  ( (k+2)(k-m  +  l,(k+m  +  l)/ (k+1  '(2k+ 1,(2 k+3)]  1/  2 

x  g(p,  k+1,  m,  A)J  -±f”  pdpj^'pr) 

x  J^(k+2,|(k-l)(k-m)(k+m)/k(2k-l)(2k+l,j  X/2  g(p,k-l,m,x) 


+  (k-1  )(  (k+2 )  (k  -  in  + 1 )  (k+m +l)/(k+l)(2k+l)(2k+3)j  1 /Z 

v  g(p,k-*l,  m,x)J  |  . 


(156) 


!  i 


* 


co  k 

>J/2  V  V  (1>kvm.0i 


sin  0afl-uo(x)  -  (2/tt)1'2  £  £  (i)K  0  (9.  *) 


k=0  m*-k 


*  |  ~2  dP  3k(pr)  ^<k-  l)(k+l)[  (k-m)(k+m)/(2k-l)(2k+l>]  l/2 

x  g(p,k-l.rn,0)  +  k(k+2)[  (k-m+l)(k+m+l)/(2k+l)(2k+3)J  1 /2 
x  g(p,  k+1,  m,  0)j  +  dP  P  jjjfer) 

x  |Wl)[<k-m)(k+m)/(2k-l)(2k+l)]  1/2g(p,  k- 1,  m,  0) 

-  (k+2)[(k-m+l)(k+m+l)/(2k+l)(2k+3)]  1/2g(p.k+l,m,0)J  |  .  (157) 


Finally, 


rv 

8in  9  <SL>  =  -<lr>‘^2£  £  (l)k[k(k+l)j'1/2Y™'°(0,*) 


k=0  m»-k 


*  |f  JQ  P  dP  Vpr)  [m8(p»  k.m.X)  +  X(k+1) 

x  f  (k-l)(k+lHk-m)(k+m)/(2k-l)(2k+l)]1/2g(p,k-l,m,X.) 

+  Xk(k(k+2)(k-m+l)(k+ni+l)/(2k+l)(2k+3)]  1/?g(p,  k+1, m, X)J 

+  J  P2dp  j^pr)  Jmg(p,  k,  m, X) 

+  X[  (k-  l)(k+l)(k-m)(k+m)/(2k-  l)(2k+l)J  k- 1,  m,  X) 

-  X[k(k+2)(k-m+l)(k+m+l)/(2k+l)(2k+3)j  1,/2g(p,  k+1,  m. X)j  j.  (158) 


lor  X  ■  ±  1,  and 


oo  k 


8in9li,ioy,'|2'/lr)l/2F  I  £  (i)k">Y™’°(0,i 


k«0  m «-k 


P  dp  jk(pr)  g(p,  k,  m,  0)  . 


(159) 
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In  Eq.  (158),  the  ambiguity  of  the  term  for  k  =  0  is  evd’uated  in  the  obvious 

way  by  setting  g  (p, -l,n,X)  =  gfa,  0,m,X)  =  0  and  then  by  multiplying  through  bv 
-l/2  ' ! 2 

the  factor  [k(k+l)J  '  (so  as  to  cancel  the  factor  k  ')  before  setting  k  =  0  . 

Thus  our  program  of  exhibiting  the  components  of  the  rotational  and  irrota- 
tional  parts  of  an  arbitrary  vector  in  a  spherical  coordinate  system  has  been 
completed.  It  is  interesting  to  note  that  while  general  surface  harmonicB  YtJ‘,n 
had  to  be  introduced  to  define  g(p,j,m,x).  the  expansions  of  the  components  of  the 
vectors  involve  the  ordinary  spherical  harmonics  Y”?'®.  The  components  of  the 
curl  01  tiie  vector  in  spherical  coordinates  have  the  same  expansion  except  inat 
g(p,  j, m, X)  is  replaced  by  Xpg(p,  j, m, X). 
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Appendix  A 

Motivation  for  tho  Introduction  ot  tho 
Eigtnfuoctiono  ol  th«  Curl  Oporatoi 


It  will  now  be  shown  how  we  came  to  use  the  vectors  Q,  (rj  >  from  considera- 
tions  of  properties  of  representations  of  the  rotation  group. 

Let  us  define  the  vector  G(j^)  as  being  the  Fourier  transform  of  the  vector 

u(x): 


G  (g_)  -  (2»)"3'f2y'dx  e'^'l  u(x)  . 
u^(x  )  *  (2ff)  “J djo  e'£-'  i  G(g_)  . 


(Al) 

(A2) 


Then 


Vtu(s)  «  i(2tr) 


-S/2 


/• 


fe  e  ^ 


[nv  G(g^)]  . 


(A3) 


Let  as  now  introduce  the  column  vector  G(g^)  whose  components  from  top  to 
bottom  are  the  x-,  y-.  and  x.~  components  respectively  of  G(p^).  It  is  readily 
shown  that  the  column  vector  corresponding  to  the  vector  jjxG(j^)  is  -i(p-  K)G(g), 
where  p  -  K  »  2  p.K.  ana  the  matrices  K.  are  given  bv  Eq.  (125).  As  is  shown  in 
Eq.  24)j  these  matrices  are  the  infinitesimal  generators  of  the  rotation  matrix. 
Let  us  now  define  the  column  vector  li(p  )  by 


* 


A  2 


G(p  j  *  VH(j) )  , 

H(p  )  ■  V  *G(£)  ,  (A4) 


where  V  is  the  matrix  oif  Eq.  (.127).  Then  the  column  vector  (p  •  K }G (p )  maps  into 
the  column  vector  <p  .3  '  ')H(p),  where  the  matrices  S:  have  the  elements  given 
by'  (106).  This  mapping  fallows  from  Eq.  (126)..  Because  of  the  way  that  the; 
row »  a. net  columns  of  the  representation*:  of  the  generators  of  the  rotation  group 
were  hi. belled,  the  top  component  of  H(jo)  will  be  labelled  H ^ (j-A,  the  second 
component,  by  H^{p  ),  and  the  hotter/:  component  by  H  ,  (p  ).  Now  In  Moses  (1970) 

and  in  others  mentioned  in  that  reference,  a  representation  is  introduced  culled 

(1) 

the  helicity  representation,  in  which  p-S  is  diagonal.  Thus  let  vis  introduce 
the  column  vector  g(p  )  whose  top,  middle,  and  hotter,/  components  are  denoted 
by  gj<p  ).  g0(p  )  and  g_.(p  ),  respectively  by 


g(p )  ■  e 


o  (1) 

~  H(p). 


H(p> 


iivj-  S 
e  — 


(1) 


«<P,> 


(A5) 


where  u  is  a  function  ctl!  grgiven  by 


p  a  -p  -a  sin  u  ,  p„  «  p  —  sin  w  ,  p„  *  p  cos  w  , 

1  vO  ■!  it)  J 

w„  «  0  ,  W  ■  |  w  |  .  (A5a) 

O 


The  relation  olu  to  the  polar  coordinates  of  p  is  given  by  Eq.  (130). 

^  / 1  \  ~~ 

The  column  vector  (p*S  '  )H{p)  maps  into  the  column  vector  g*(p>  whose 
components  are  ^Pfex<E>'  Ttlus«  combining  all  tile  transformations,  u_  (x  )  maps  into 
the  column  vector  g(js )  in  such  a  way  that  V  x  u  (x  )  maps  into  the  column  vector 
g'(p )  .  Equations  (13)  are  an  explicit  form  for  the  mapping  and  Eqs.  (12)  and 
(129)  constitute  a  statement  of  the  mapping  of  the  curl  operator. 
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